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摘 要 : 博 斯 腾 湖 径流 入 湖 、 出 湖 口 过 近 的 水 文 现状 对 有 机 物 浓度 及 其 空间 分 布 形成 制约 。 采 用 MIKE 1 软件 构 


建 水 动力 水 质 二 维 数值 模型 ,分 析 工 程 方案 对 博 斯 腾 湖 有 机 物 污染 的 改善 效果 。 根 据 模拟 结果 以 及 配对 样本 7 检 


验 分 析 : 在 现状 基础 上 实施 扬 水 站 东 移 或 引 开 都 河水 3 x10 m° 


“a "进入 黄 水 沟 , 均 可 显著 改善 博 斯 腾 湖 有 机 物 污 


染 , 年 均 化 学 需 氧 量 (COD) 浓 度 分 别 降低 13.4% 和 8.91% ;两 方案 同时 实施 效果 更 为 显著 ,年 均 COD 浓度 降低 


20.2% ;通过 工程 实施 改善 博 斯 腾 湖 有 机 物 污染 切实 可 行 。 但 在 流域 尺度 上 ,工程 方案 仅 改 变 了 有 机 物 平衡 而 未 


关键 词 : 博 斯 


专 斯 腾 湖 是 我 国 最 大 内 陆 湖泊 ,在 区 域 社会 经 
济 发 展 . 生 态 平 衡 维持 ,生物 多 样 性 保持 等 方面 发 挥 
着 重要 作用 "“”。 随 着 矿 化 度 上 升 有 机 物 污染 .水 
资源 短缺 等 问题 凸显 , 博 斯 腾 湖 生态 环境 问题 受 
到 学 者 关注 ,包括 围绕 水 质 水 量 ”“ k E 
物 “ ”水 生生 物 "“ DR O ”的 基础 研究 , 面 
向 湖泊 需 水 量 "” AE AS KB aK BEBE) EE 
理 研究 以 及 针对 矿 化 度 分 布 “ \ 塔 里 木 河 输 水 
等 实际 问题 的 工程 探讨 。 
章 斯 腾 湖 核心 环境 问题 在 于 有 机 物 污染 ,突出 
表现 为 大 部 分 湖区 水 体 化 学 需 氧 量 ( COD , chemical 
oxygen demand) 浓 度 高 于 地 表 水 亚 类 标准 限 值 ,以 
及 空间 差异 显著 “ 。 有 机 物 对 湖泊 自净 能 力 反映 
较 弱 ” ,浓度 降低 更 多 依靠 稀释 交换 作用 ,而 博 斯 
腾 湖 径流 入 湖 、 出 湖 口 位 置 较 近 ,径流 对 湖水 稀释 作 
用 更 多 体现 在 西南 局 部 湖区 。 在 现状 水 文 条 件 下 ， 
大 部 分 湖区 有 机 物 无 法 得 到 稀释 ,也 难以 通过 水 体 
交换 而 排出 ,阻碍 博 斯 腾 湖 全 面 水 质 达 标 。 因 此 , 通 
过 工程 优化 径流 入 湖 及 出 湖 口 ,继而 促进 湖水 稀释 
交换 成 为 改善 博 斯 腾 湖 有 机 物 污 染 的 潜在 方案 。 

本 文采 用 MIKE21 软件 构建 水 动力 水 质 二 维 数 
值 模型 ,模拟 不 同 工 程 方案 对 于 博 斯 腾 湖 有 机 物 污 
染 的 改善 效果 ,以 期 为 博 斯 腾 湖 水 环境 保护 提供 理 
论 依据 和 工程 参考 。 
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溺 湖 有 机 物 污染 改善 源 于 污染 物 转移 至 下 游 孔 淮河 。 
BW: 有 机 物 污染 ; MIKE21; COD; 水 文 条 件 


1 研究 区 域 概况 


专 斯 腾 湖 位 于 中 国 天 山南 节 在 普 伟 地 东南 部 低 
洼 处 ,地 理 位 置 介 于 86°19’ ~ 87°28’E, 41°46’ ~ 
42°08'N ,是 开 都 河 等 水 系 的 尾 间 以 及 孔 和 淮河 源头 
(图 1)。 博 斯 腾 湖 由 公路 分 隔 为 大 小 湖区 :大 湖区 
在 水 位 1 047. 00 m 时 ,面积 1 064. 1 km’, 容积 
73.03 x 10° mi ,平均 水 深 约 7.5 mm, 最 深 处 约 16 m; 
小 湖区 为 诸多 浅 小 湖泊 和 芦苇 沼泽 湿地 。 

得 斯 腾 湖 补给 河流 历史 上 有 开 都 河 、 黄 水 沟 清 
水 河 、 曲 囊 沟 、 乌 拉 斯 台 河 ,现状 仅 余 开 都 河东 支 。 
湖滨 的 黄 水 沟 湿地 承担 着 污水 深度 净化 功能 ,部 分 
污水 通过 渗 漏 以 及 汛期 行 洪 进入 大 湖区 。 博 斯 腾 湖 
唯一 出 水 口 是 位 于 大 湖区 西南 角 的 扬 水 站 ,水 体 由 
东 \ 西 两 个 泵 站 提升 进入 孔 稚 河 。 


2 MIKE21 模型 


2.1 基本 方程 

MIKE21 是 丹麦 水 利 研 究 所 开发 的 系列 水 动力 
学 软件 之 一 ,在 水 环境 模拟 中 得 到 诸多 应 用 。 
专 斯 腾 湖 为 浅水 湖泊 ,水 平 尺 度 远大 于 垂直 尺度 ,用 
MIKE21 HD 模块 构建 垂 向 平均 的 二 维 水 动力 模型 ， 
控制 方程 如 式 (1) 、 式 (2)、 式 (3); 采 用 MIKE21 
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图 1 博 斯 腾 湖 水 系 图 
Fig.1 The water system of Bosten Lake 


Transport 模块 构建 水 质 模型 ,基本 方程 如 式 (4) 。 
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式 中 :i 为 时 间 ;u 及 vw 为 流速 ;g 为 重力 加 速度 ;m 为 
水 位 ;h 为 水 深 ;f 为 科 氏 力 系 数 ;p 为 密度 ,S 为 源 
项 ;rs rw n Tin Ty N y 方向 的 表面 风 应 力 和 底部 
切 应 力 。 
OO «= a) 
ot Ox oy i 
式 中 :Ff 为 水 平 扩散 项 ;C 为 垂 向 平均 浓度 ;$, 为 源 
CTH sk, 为 线性 衰减 速率 ;C, 为 污染 物 浓度 。 模 型 
计算 范围 为 博 斯 腾 湖 大 湖区 ,采用 三 角形 不 规则 网 
格 ( 图 2) ,三 角形 边 长 约 为 50 m, 
2.2 模型 验证 
从 新 疆 巴 音 郭 楞 蒙古 自治 州 环 境 保 护 局 收集 
2012 一 2016 年 博 斯 腾 湖 17 个 国家 地 表 水 环境 质量 
监测 点 水 质数 据 、 开 都 河东 支 水 质数 据 及 污染 物 数 


Bathymotry/m 
= Above~-0.6 =-1.6~-0.8 = -2.4~-1.6 = -3.2~-2.4 
= —4,0~-3.2 =-4,.8~-40 =-5.6~-48 = -6.4~-5.8 
= -7.2~-6.4 = -8.0~-7.2 =-8.8~-8.0 = -9.6~-8.8 
= -10.4~-9.6 =-11.2~-10.4 =-12.0~-11.2 = Below~12.0 


= Undorned Vaue 


2 模型 计算 网 格 
Fig.2 The computational grid of MIKE21 


据 , 从 塔里木 河流 域 管理 局 收集 开 都 河 2013—2016 
年 东 文 径流 量 \ 扬 水 站 输 水 量 . 博 斯 腾 湖 水 位 等 月 均 
水 文 数据 。 鉴 于 博 斯 腾 湖 湖 流 受 风 场 影 响 较 大 的 特 
HE ,从 气象 局 收集 2013—2016 年 降雨 .蒸发 、 风 
速 .风向 等 月 均 气 象 数据 作为 模型 基础 输入 参数 。 

以 2013 一 2015 年 数据 进行 模型 率 定 , 以 2016 
年 数据 进行 结果 验证 。 模 型 初始 条 件 为 2013 年 1 
月 1 日 ,时 间 步 长 为 5 min ,初始 水 位 为 实测 水 位 , 初 
始 COD 浓度 以 2012 年 10 月 实测 浓度 代替 。 

模型 模拟 值 与 实测 值 对 比如 图 3 所 示 , 采 用 过 
程 线 拟 合 中 常用 的 均 方 根 误差 (RMSE) 纳什 效率 
系数 (NSE ) 两 个 函数 进行 验证 “1 。 RMSE 计 算 
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Fig.3 Verification of water level and COD 


函数 如 式 (5) ,表征 模拟 值 与 实测 值 之 间 的 偏差 。 
NSE 计算 函数 如 式 (6) ,表征 两 组 数据 序列 的 模拟 
精度 :接近 1 则 表示 模拟 质量 越 好 、 模 型 可 信和 度 越 
高 ;接近 0 则 表示 模型 总 体 结果 可 信 , 但 过 程 模拟 误 
差 较 大 ;小 于 0 则 表明 模拟 结果 不 可 信 。 


12 i i\2 

RMSE =, /--3 (Qi - 0) (5) 
O 

NSE =1 -后 一 一 一 (6) 
ECO - 0)” 


RP: Q 为 实测 值 ;0, 为 模拟 值 ; Oo 为 实测 值 的 平 
均值 。 

验证 表明 ,水 位 模拟 结果 良好 (RMSE = 0. 056, 
NSE = 0. 97) , COD 模拟 结果 整体 可 信 (RMSE = 
1.77,NSE =0.36) 。 


3 ”结果 与 讨论 
3.1 工程 方案 

针对 博 斯 腾 湖 径流 入 湖 、 出 湖 口 距离 较 近 的 状 
况 ,结合 区 域 实际 情况 和 工程 可 行 性 ,设计 3 类 工程 
方案 ( 表 1) :径流 入 湖 优化 、 径 流出 湖 优化 .入 湖 出 


湖 同 时 优化 。 

径流 出 湖 优 化 从 增强 开 都 河 对 湖水 稀释 作用 、 
提高 出 湖水 体 污染 物 总 量 角 度 出 发 ,将 现状 扬 水 站 
向 东部 迁移 ,根据 湖岸 特征 设计 3 个 迁移 位 置 :1# 为 
南岸 最 北部 ;2# 为 二 道 海 心 山 附近 ;3# 为 三 道 海 心 山 
附近 。 博 斯 腾 湖 南部 现状 为 荒漠 区 域 , 扬 水 站 迁移 
及 新 建 输 水 设施 均 具 有 较 高 可 行 性 。 

径流 入 湖 优 化 从 恢复 博 斯 腾 湖 自然 水 文 条 件 出 
发 ,将 开 都 河东 支部 分 径流 引入 黄 水 沟 ( 多 年 平均 
径流 量 为 2.91 x 10° m*) ,设计 引水 径流 量 分 别 为 
1x 10° 2 x10° 3x10 m+ a 的 3 种 方案 。 开 都 
河 与 黄 水 沟 均 位 于 珊 普 盆地 ,相互 距离 较 近 且 地 势 平 
坦 ,引水 工程 具有 可 行 性 。 此 外 ,引水 工程 可 作为 黄 
水 沟 汛期 洪水 入 湖 通道 ,从 而 减少 湿地 污水 入 湖 量 。 

入 湖 出 湖 同时 优化 则 为 上 述 两 类 方案 中 最 优 方 
案 的 组 合 ( 方 案 2+ 方 案 6)。 
3.2 模拟 结果 

各 方案 2016 年 模拟 结果 统计 分 析 如 表 2。 在 
收集 2012—2016 年 的 博 斯 腾 湖 水 质数 据 中 ,最 后 一 
批 数据 为 2016 年 10 月 15 日 ,各 方案 模拟 的 该 日 
COD 浓度 场 如 图 4。 


表 1 工程 方案 设计 
Tab.1 Design of projects 


工程 方案 开 都 河东 支 年 径流 量 黄 水 沟 年 径流 量 /10? m? 扬 水 站 位 置 黄 水 沟 湿 地 污染 物 输 入 
现状 现状 值 0 西南 角 模型 率 定 值 
径流 出 湖 优化 方案 1 现状 值 0 东南 角 ]# 模型 率 定 值 
方案 2 现状 值 0 东南 角 2# 模型 率 定 值 
方案 3 现状 值 0 东南 角 3# 模型 率 定 值 
径流 入 湖 优化 方案 4 现状 值 减少 1 x 108 m? 1 西南 角 减少 至 率 定 值 1/5 
方案 5 现状 值 减少 2 x 108 m3 2 西南 角 减少 至 率 定 值 1/5 
方案 6 现状 值 减少 3 x108 m3 3 西南 角 减少 至 率 定 值 L/S 
同时 优化 方案 7 现状 值 减少 3 x108 m3 3 东南 角 2# 减少 至 率 定 值 1/5 


ChinaXiv 合 作 期 刊 


胡 春 明 等 : 博 斯 腾 湖 有 机 物 污染 改善 方案 研究 


表 2 各 方案 2016 年 化 学 需 氧 量 ( COD ) 浓度 模拟 结果 统计 
Tab.2 Simulation of COD in 2016 


COD 浓度 / (mg + L-!) 
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工程 方案 
5 月 6 月 7 月 8 月 9 月 10 月 年 均 
现状 22.26 +1.29 22.16 +1.78 21.85 +1.47 22.18 +1.64 21.96 +1. 16 21.89 £0.87 22.05 £1.19 
方案 1 19.60 +1.52 19.58 +1. 86 19.20 +1.75 19.05 +2.27 18.70 +1.61 18.77 £1.04 19.15 +1.46 
方案 2 19.54 +1.55 19.52 £1.87 19.13 +1.76 18.99 +2. 26 18.67 +1.61 18.77 £1.03 19.11 +1.45 
方案 3 19.63 +1.54 19.60 +1.87 19.20 +1.77 19.05 +2.27 18.72 +1.61 18.83 +1.04 19.17 £1.45 


方案 4 20.84 +1.06 20.49 +1.28 20.11 +1.27 20.28 +1.48 20.20 +0.92 20.00 +0. 92 20.32 +0. 97 
方案 5 20.73 +1.08 20.37 +1.30 19.99 +1.35 20.17 +1.55 20.11 +0.91 19.91 +0.91 20.21 +1.00 
方案 6 20.60 +1.10 20.23 +1.32 19.84 +1.45 20.04 +1.63 20.00 +0.92 19.79 +0.91 20.09 +1.05 
方案 7 18.27 +1.24 17.96 +1.29 17.52 +1.45 17.35 41.75 17.24 41.13 17.19 +0.92 17.59 +1.07 


注 :数据 为 “均值 + 标准 差 ”。 
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于 时 


径流 出 湖 优化 的 3 种 方案 结果 相差 不 大 ,分别 
可 实现 博 斯 腾 湖 年 均 COD 浓度 下 降 13. 2% 、 
13.4% 和 13.0% 。 径 流入 湖 优化 的 3 种 方案 分 别 
可 降低 年 均 COD 浓度 7.84% .8.33% 和 8.91% , 黄 
水 沟 年 径流 量 越 大 ,湖区 水 质 改善 效果 愈 好 。 而 径 
流入 湖 出 湖 同时 优化 则 可 实现 博 斯 腾 湖 年 均 COD 
浓度 下 降 20.2% 。 

基于 模拟 结果 的 标准 差分 析 ,方案 1 方案 2 N 
K 3 模拟 结果 的 标准 差 较 现状 略 有 上 升 ,各 监测 点 
COD 浓度 的 离散 程度 更 大 ,表明 博 斯 腾 湖 COD Ye 
度 整 体 空间 差异 性 增 大 ,其 原因 在 于 出 湖 口 优化 至 
东南 角 , 仅 改善 了 南部 区 域 COD 浓度 , 却 进 一 步 加 
大 了 与 北部 区 域 的 差异 。 其 他 方案 模拟 结果 的 标准 
差 总 体 下 降 , COD 浓度 空间 差异 性 减 小 ,关键 原因 
在 于 径流 从 黄 水 沟 入 湖 降低 了 现状 污染 最 为 严重 的 
西北 区 域 COD 浓度 。 
3.3 方案 检验 

在 各 方案 模拟 结果 中 ,方案 2 是 径流 出 湖 优 化 
的 最 佳 方案 ,方案 6 是 径流 入 湖 优化 的 最 佳 方案 ,而 
方案 7 则 为 两 者 的 组 合 。 采 用 “配对 样本 7 了 检验 ” 
分 析 方 案 实施 对 博 斯 腾 湖 水 质 改善 效果 ,根据 检验 
结果 ( 表 3) ,各 配对 样本 检验 结果 对 应 的 P<0.05， 
表明 在 现状 基础 上 实施 方案 2 或 方案 6 均 可 显著 改 
善 博 斯 腾 湖 COD 浓度 ,而 最 终 实施 方案 7 则 可 实现 
COD 浓度 进一步 显著 改善 。 


RI 各 方案 配对 样本 了 检验 结果 
Tab.3 Results of paired samples 7-test 


现状 - 现状 - 方案 2 - 方案 6 - 

方案 2 KRE KRIT VRI 

7 统计 值 57.363 11. 684 9. 112 56. 204 
相伴 概率 (P) 0.000 0. 000 0. 000 0. 000 


3.4 讨论 

湖泊 污染 物 浓度 上 升 根本 原因 在 于 总 量 不 平 
衡 , 入 湖 量 高 于 出 湖 量 带 来 污染 物 累 积 。 针 对 博 斯 
腾 湖 矿 化 度 问 题 , 有 研究 提出 优化 设想 : 韩 龙 喜 
等 "5 模拟 扬 水 站 不 同位 置 对 矿 化 度 空间 分 布 产 生 
影响 ; 吴 敬 禄 等 分析 了 扬 水 站 建设 运行 使 博 斯 腾 
湖 总 盐 量 基本 保持 在 9.60 x 10° + 的 平衡 点 ,而 扬 水 
站 东 迁 可 将 平衡 点 降低 在 8.50 x 10° to 

湖泊 污染 治理 通常 采用 污染 源 系统 治理 、 清 水 
产 流 机 制 修复 、 湖 泊 水 体 生 境 改善 .流域 管理 等 措 


fa 


wit) 。 随 着 博 斯 腾 湖 流域 污染 源 治理 及 生态 保护 
工程 不 断 实 施 ,湖泊 水 质 持 续 改善 ,年 均 COD 浓度 
由 2013 年 28. 68 mg * L-! 降 至 2016 年 21.50 mg - 
Lo! COD 浓度 也 逐渐 达到 一 个 平衡 点 。 若 要 继续 
改善 湖区 水 质 ,提高 污染 源 治理 及 生态 保护 工程 效 
果 , 需 通过 工程 手段 改变 博 斯 腾 湖 水 文 条 件 .打破 有 
机 污染 物 平衡 。 

开 都 河水 体 COD 浓度 普遍 在 10 mg ' 工 ”以 
内 , 汇 人 后 在 西南 区 域 与 湖水 混合 ,再 由 扬 水 站 排 和 
下 游 孔 汰 河 ,出 湖水 体 携 带 的 有 机 物 总 量 较 少 。 本 
研究 设 定 的 各 工程 方案 本 质 均 为 增加 开 都 河 径流 停 
留 时 间 ,提高 径流 与 湖水 稀释 混合 程度 ,打破 博 斯 腾 
湖 有 机 物 现状 平衡 。 博 斯 腾 湖 水 体 的 5 日 生化 需 氧 
量 (BOD5 ) 浓度 低 ( 约 2 mg . Lo!) 浮游 生物 较 
少 ,有 机 物 生物 降解 较 弱 ,稀释 交换 作用 对 有 机 
物 浓度 影响 尤为 显著 ,使 各 方案 取得 较 好 改善 效果 。 
但 从 流域 尺度 考虑 ,各 方案 并 未 减少 有 机 物 总 量 , 仅 
改变 有 机 物 时 空 分 布 , 博 斯 腾 湖 有 机 物 总 量 降低 , 则 
意味 着 下 游 孔 省 河 有 机 物 总 量 升 高 ,生态 环境 将 受 


到 一 定 影 响 。 
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4 结论 


径流 入 湖 .出 湖 口 位 置 较 近 的 水 文 现 状 对 博 斯 
腾 湖 COD 浓度 及 其 空间 分 布 形 成 制约 。 通 过 优化 
径流 入 湖 、 出 湖 可 改变 博 斯 腾 湖 有 机 物 平衡 ,在 不 改 
变 有 机 物 总 量 的 情况 下 降低 湖区 COD 平衡 浓度 :将 
扬 水 站 东 移 可 实现 全 湖 年 均 COD 浓度 较 现状 降低 
13.4% , 引 开 都 河水 3 x 10° m”. a ! 进 入 黄 水 沟 可 
降低 年 均 COD 浓度 8. 91% ,两 者 同时 实施 可 降低 
年 均 COD 浓度 20. 2% ;后 两 种 方案 还 对 湖区 COD 
浓度 空间 差异 性 也 有 一 定 改善 作用 。 
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Study on the project of reducing organic pollutants in the Bosten Lake 


HU Chun-ming', Narengerile’, MA Jin-feng', QING Jin', YOU Li’ 
(1. Research Center for Eco-Environmental Sciences , Chinese Academy of Sciences , Beijing 100085 , China; 
2. Institute of Lake Bosten ,Bayingolin Mongolia Autonomous Prefecture of Xinjiang ,Korle 841000 , Xinjiang , China) 


Abstract; The contiguous river estuary and outfall influence the discharge of organic pollutants, which limit the 
concentration and spatial distribution of the organic pollutants in the Bosten Lake. Various projects should be con- 
ducted to optimize the current hydrological conditions. The MIKE21 software was applied to establish a two-dimen- 
sional numerical model and analyze the effects of different intervention projects. According to the simulation results 
and the paired sample T-test results ,the annual chemical oxygen demand concentration in the Bosten Lake could be 
reduced by 13.4% by moving the outfall from the southwest to the southeast 8.91% by diverting 3 x 10% m° + a7! 
of the runoff from the Kaidu River to the Huangshuigou River, and 20.2% by implementing both these projects. The 
results denote that it is feasible to reduce the organic pollutants in the Bosten Lake by optimizing the hydrological 
conditions. However ,these projects only changed the proportion of the organic pollutants in the Bosten Lake, they 
did not reduce the quantity of organic pollutants in the basin. The organic pollutants were diverted into the Kongque 


River downstream. 


Key words; Bosten Lake; organic pollutants; MIKE21; chemical oxygen demand; hydrological condition 


